Steady-state superradiance with alkaline earth atoms 
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Earth-alkaline-like atoms with ultra-narrow transitions open the door to a new regime of cavity 
quantum electrodynamics. That regime is characterized by a critical photon number that is many 
orders of magnitude smaller than what can be achieved in conventional systems. We show that it is 
possible to achieve superradiance in steady state with such systems. We discuss the basic underlying 
mechanisms as well as the key experimental requirements. 

PACS numbers: 42.50.Fx, 37.30.+i, 42.50.Pq, 42.50.Ct, 42.55.Ah, 42.50.Lc 



Superradiance, first introduced by Dicke over 50 years 
ago p| , is one of the pillars of cavity quantum electrody- 
namics (CQED). Superradiance occurs due to the con- 
structive interference of the probability amplitudes for 
spontaneous decay of several atoms. Due to its general- 
ity and conceptual simplicity it is a paradigm system for 
collective behavior. Usually superradiance is transient; 
atoms initially prepared in the excited state relax to the 
ground state rapidly and the collective emission termi- 
nates. To date superradiance has not been achieved in 
a continuous fashion. The goal of this letter is to show 
that steady-state superradiance can be achieved with ul- 
tracold alkaline earth atoms in high finesse cavities. Such 
systems are experimentally available in the form of opti- 
cal lattice clocks @, Q . Recently Bose Einstein conden- 
sates of such atoms have also become available with Ca 
and Sr atoms 

Atoms with a two-electron level structure possess nar- 
row inter-combination lines that, due to selection rules, 
are dipole-forbidden, and typically have lifetimes many 
orders of magnitude longer than dipolc-allowed transi- 
tions. The long lived excited state is essential for super- 
radiance in steady-state because it allows the buildup of 
population inversion even as the population of the ex- 
cited state is drained by the collective decay. For dipolc- 
allowed transitions, on the other hand, supcrradiant de- 
cay is so rapid that it would exhaust the supply of excited 
state atoms before it could be replenished by repumping, 
and consequently the superradiant emission must cease. 

Besides being of fundamental importance, steady state 
superradiant systems are also interesting because of their 
potential applications. The most immediate application 
is the possibility to build active optical clocks where the 
light serving as a frequency standard is derived directly 
from the atoms 0, H[ ■ Such systems have the potential 
to improve the stability of the best clocks by about two 
orders of magnitude. Another area of application is to 
strongly-correlated physics. The atoms in this system 
evolve into exotic many-particle states due to their col- 
lective interaction with the light field. These states could 
be of interest for quantum information purposes as well 
as for the exploration and study of many-particle phases 
of condensed matter. 

The model that we consider is depicted in Fig.QJa). N 
two level atoms with excited state |e) and ground state 



FIG. 1: (Color online) (a) Schematic of a coupled atom-cavity 
system leading to steady state superradiance. (b) The param- 
eter space of CQED spanned by the critical atom number no 
and the critical photon number mo. The collective behavior of 
the coupled atom-cavity system undergoes a cross-over from 
stimulated emission dominated (laser-like) at large mo to col- 
lective spontaneous emission dominated (superradiance-like) 
for small mo. 



\g) decay with rate r c = CT through the mode of a cavity, 
where C = g 2 /(Tn) is the single atom coopcrativity pa- 
rameter. The transition from \e) to \g) is assumed to be 
a narrow inter-combination line. The atomic free space 
spontaneous emission rate is T, the single photon Rabi 
frequency is g, and the cavity decay rate is k. For sim- 
plicity we assume that all atoms couple identically to the 
cavity mode. The rate T c is assumed to be much larger 
than the rates for non-collective decay processes and de- 
phasing. That approximation requires that NC 3> 1. At 
the same time the atoms are being non- collectively re- 
pumped to the excited state with an effective rate w. 
This could be achieved by resonantly driving a transi- 
tion to a third state with Rabi frequency fl that decays 
rapidly with rate 7 to |e). 
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There are two key experimental requirements for re- 
alizing the physics discussed here. The most impor- 
tant requirement is that the collective decay be domi- 
nant over all other decay and decoherence processes, i.e. 
NCT ^> r, I/T2, where T2 is the spin dephasing time. 
The second requirement is that one must be able to re- 
populate the atomic excited state at a rate equal to the 
collective decay rate, i.e. it must be possible to achieve 
w ~ NCT without significant atomic losses. 

The model under consideration is similar to typical 
many-atom CQED systems that are studied theoreti- 
cally and experimentally, except that the atomic dipole 
moment is many orders of magnitude smaller. Never- 
theless, the small dipole moment can lead to profound 
consequences. These can be characterized by the criti- 
cal atom number no = (nT)/(g 2 ) and the critical pho- 
ton number itiq = (7 2 )/(g 2 ) that tell us to what degree 
quantum effects arc important: no < 1 means that a 
single atom can substantially affect the cavity field and 
tuq < 1 means that the electrical field corresponding to a 
single photon in the cavity can saturate the atomic tran- 
sition. For ultra-narrow inter-combination lines in earth- 
alkalinc-like atoms, critical photon numbers as small as 
mo ~ 10" 12 can realistically be achieved. These sys- 
tems reside in an exotic region of parameter space that 
has previously been inaccessible. For example for cavity 
QED systems with alkali atoms the records are in the 
mo, n a ~ 0.01 — 0.001 range. It is impossible to sub- 
stantially improve upon these values as can be seen by 
rewriting the critical photon number as no = 2nA/(Fa) 
and Too = 47r 2 V r c ff/(<5Ao) where A is the cross section 
of the cavity mode, a = 3A 2 /(27r) is the resonant cross 
section of the atoms, Ao is the wavelength of the reso- 
nant light, and F and Q are the finesse of the cavity 
and the quality factor of the atomic transition, respec- 
tively. Since the dimensionless V c s/(Xq) must be at least 
unity on fundamental grounds, and is typically orders of 
magnitude larger, the atomic resonance Q must be ex- 
traordinarily high to reach mo ~ 10 -12 , and the values 
of Q that are possible for the optically allowed dipole 
transitions in alkali atoms and similar systems are insuf- 
ficient. Incidentally, masers, which operate in the mi- 
crowave domain, have critical photon numbers similar to 
the narrow linewidth atoms, but are classical due to their 
rather large critical atom number, and have photon en- 
ergies that are many orders of magnitude smaller than is 
characteristic in the optical domain. 

The evolution of the system shown in Fig. [T] is given 
by the master equation, 

^ = - y (jJ-P + pj+J- - 2J-/5J+) (1) 

jv 

- f £ + - 2*?/»*- W ) • 

3=1 

Here, a^P = (a+ )^ = |fif)(e| is the lowering operator for 

atom j and J_ = (J+)t = J2f=i 1S tne collective de- 
cay operator brought about by the coupling of the atoms 
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FIG. 2: (Color online) Emission rate as a function of repump 
rate for 10 atoms. The orange dots are the expectation values 
extracted from the Monte-Carlo wave function simulations. 
The blue line shows the semi-classical result obtained from 
Eqs.(|3][5]), and the purple dashed line is the emission rate for 
uncorrelated particles, N e T c . 

to the rapidly decaying cavity field. In deriving this mas- 
ter equation we have assumed that the cavity field decays 
so rapidly that the mean photon number is much smaller 
than unity, and thus the field has been adiabatically elim- 
inated. This approximation is valid for atoms with an 
extremely weak dipole moment such as those that we are 
interested in. The collective decay part occurs in the 
superradiance master equation first introduced by Boni- 
facio et. al. [H; see also The repumping can be 
thought of as spontaneous "absorption" from the ground 
state to the excited state. 

We unravel the master equation into Monte Carlo 
wavefunction trajectories [ll| . and from the trajec- 
tories \ip(t)) we extract expectation values (0(t)) = 
(ip(t)\0\ip(t)) of system observables O. Steady state ex- 
pectation values (0)ss are then obtained by calculating 
time averages, (O)ss = ^ / t * 1+T di (0(f)) , where t\ is 
chosen large enough to allow the system to settle to 
steady state, and T is chosen long enough for statisti- 
cal errors to be controlled. We have empirically checked 
that the steady state does not depend on the choice of 
initial conditions. 

In Fig. [5] we show that Eq. |T]) leads to sustained su- 
perradiance by investigating the mean photon emission 
rate 

I = T c (J+J-)ss, (2) 

together with the emission rate A^ e r c that one would ex- 
pect for uncorrelated atoms, with N e the population of 
the excited state. Three qualitatively different regimes 
can be distinguished: strong pumping with w > NT C , 
intermediate pumping with T c < w < NT C , and weak 
pumping w < T c . For very strong pumping, w 3> iVr c 
the emission rate approaches the maximum possible emis- 
sion rate for uncorrelated atoms, -/Vr c . As the pump 
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FIG. 3: (Color online) Transition rates between J and M 
eigenspaces (a-c) and populations Pm,.j of these subspaces 
(d-f) for w = 0.1 T c (a, d), to = 2.0 T c (b, e), w = 10.0 F c (c, 
f), and N = 4 atoms. Red arrows indicate decay dominated 
transitions and blue arrows indicate repump dominated tran- 
sitions. The transition rates are calculated by averaging over 
degenerate initial states and summing over final states. In 
(d-f) darker shades indicate larger population. 



rate decreases, the emission rate increases beyond NT C . 
This implies the presence of correlations between dif- 
ferent atoms, resulting in a collective enhancement of 
the emission rate. The behavior of the system in the 
weak pumping regime, w < T c is surprising and counter- 
intuitive in two ways. First, nearly half the atoms re- 
main in the excited state even as w/T c — > 0. Second, 
the emission rate is greatly suppressed below the value 
for uncorrelated atoms, i.e. the atoms are subradiant. 
Subradiance, like superradiance, implies the presence of 
correlations between different atoms. 

The qualitative behavior of the emission rate can be 
understood if we look at the dynamics of the system in 
the collective basis | J, M , £) where J characterizes the 
total angular momentum, the magnetic quantum num- 
ber M characterizes the inversion, and the multiplicity 
quantum number £ enumerates the degenerate manifold 
at given J and M. Figures [3ja-c) show the subspaces 
corresponding to well defined J and M along with the 
net transition rates, for N — 4 atoms. Panels (d-f) show 
the steady state probabilities, Pm,j = (Pm,j) SS where 
Pm,j = Y^e \ Jt £)(J,M,£\, for the system to be in the 
various subspaces. The right hand column (panel c and f) 
of Fig. [3] shows the strong pumping limit where rcpump- 
ing dominates the collective decay and consequently a 
large fraction of the atoms accumulate in the excited 
state. 

In the intermediate pumping regime, w < NT C , shown 



in the central column (panel b and e), there is a non- 
trivial competition between collective decay and repump- 
ing. It is apparent that for states where J ~ 0(N/2) and 
is therefore close to its maximum possible value, the col- 
lective decay dominates, whereas for states with small 
J the non-collective repumping dominates. This leads 
to a cycle in which the decay and repumping balance 
with an enhanced average value for (J + J_), i.e. an en- 
hanced emission rate. Note that it is important that the 
repumping is non-collective for superradiance to occur 
in steady state. If the repumping also preserves J, as 
is the situation considered in Refs. [H, [TH, [l4j|, there is 
no mechanism to balance the non-collective decay and 
decoherence processes that tend to decrease J and are 
unavoidable in experiment. 

In the weak pumping regime, shown in the left hand 
column (panel a and d), the collective emission drives 
the system into states with M = — J that cannot decay 
by means of J_. For M < the repumping predomi- 
nantly drives transitions J — » J— 1 while transitions from 
J — >■ J+ 1 are rare. Thus the system evolves along a dy- 
namical pathway with smaller and smaller J, eventually 
reaching J ~ 1 where it becomes trapped. This steady 
state has an almost equal number of atoms in the ground 
and excited state, and yet has a greatly suppressed emis- 
sion rate. This is due to subradiance arising from strong 
atom-atom correlations. The requirement for this subra- 
diance to occur is C 3> 1, which is much more stringent 
than the condition for superradiance. 

The Monte Carlo simulations provide a clear picture 
of the dynamical interference, and are complete in that 
they allow us to include a full description of the Hilbcrt 
space evolution. Due to the exponential scaling of the 
dimension of the Hilbcrt space of the system with the 
number of particles we are however limited to relatively 
small numbers of particles (of order 20). To consider 
mesoscopic particle numbers, we use a semi-classical ap- 
proximation that consists of keeping only pair correla- 
tions. This approximation is validated by the excellent 
agreement with the Monte-Carlo results for small particle 
numbers. Mathematically, the semi-classical approxima- 
tion is implemented by expanding expectation values of 
the system operators in terms of cumu- 

lants (. . .) c , where = |e)(e| — |.g}(.9|- We use that all 
expectation values are symmetrical with respect to parti- 
cle exchange, e.g. (a^a^) = (d+&_) for all i ^ j. All 
non-zero cumulants up to second order can be expressed 
in terms of (a { z 1} ) c = <a<V 2) ) c = <«7<V 2) >, and 

{oz Oz )c = {oz Oz ) — (o'i 1 '} 2 and their equations of 
motion are 

^^ = -( W + F c )(a«) c (3) 

-2r c (iv-i)<<7;V 2) ) c , 



4 



d(a^a^) ( 
dt 



-(™ + r c )<4V 2) ) c 
r, 



(4) 



d((Tz '(Tz ) 

dt 



+ y((^ 1) ^ ) ) c + (^ 1) ) c ) 



-2(w + T c )(a^ai% 



(5) 



By dropping the small third order cumulant 
(ai 1 ^ ) c we obtain a closed set of equations. 

The steady state is found by setting the time derivatives 
to zero. The resulting cubic equations can be solved 
exactly and these solutions are the basis of the analytical 
curve in Fig. @. The expressions one obtains are 
however very complicated. Simple expressions may be 
obtained in the regime r c <C w ~ NT C in which collective 
emission occurs. By introducing the rescaled operators 
j z = J z /N and j ± = J±/N where J z = 1/2^ , 
we find that to leading order in 1/N the only non-zero 
expectation values are (j z } c and (j+j-) c and they evolve 
according to, 



d(j z ) c 
dt 

d(j+j-)c 
dt 



-w{(j z ) c -\/2) - AT c (j + j_) c (6) 
-w(j+J-)c + 2NT c (j z ) c (J + j_) c . (7) 



The steady state solutions are 



and 



&>ss = 1^ 



w < NT C 
w > AT r 



(j + 5-)ss= 2OTo 




1 - 



wr, 



w < NT C 
w > NT C 



(8) 



(9) 



The fact that (j+j-)ss is of order unity for w < NT C 
indicates the presence of strong correlations between 
atoms. From Eq. ([9]) we can extract the maximum in- 
tensity. 



obtained at w — NT c /2. The scaling of that intensity 
with N 2 underlines the collective nature of the light emis- 
sion. Remarkably, I ma x is only a factor 1/2 smaller than 
the maximum possible emission rate r c ( J, 0| J+ J- 1 J, 0). 

Not only is this source of light bright, however, it also 
has a long coherence time — a property shared with lasers. 
In lasers, the coherence time of the light field is much 
longer than the cavity-ring-down time because of the fact 
that the field is macroscopically occupied [HI]. To find 
the coherence time of the atoms in the supcrradiant sys- 
tem considered in this paper, we study the two-time cor- 
relation function of the atomic dipole (a+\t + r)a_ (t)) . 
Using the quantum regression theorem we find the equa- 
tions of motion 



dr 



w 



(N-2)(aP)ss) 



x(^W^(0)>, (11) 



(10) 



where we have factorized (a^ {t)& z 2) (t)^ (0)) w 

(<7+'')ss(o'2 2 ' > (i")<7?' ) (C))} and we have assumed that the 

system is in steady state. Inserting (ai^ss = 2(j z )ss 
from above we find that the spins of different atoms re- 
main coherent with each other for t co h = N/(NT C + 2uu). 
In the supcrradiant regime where w ~ NT C the coherence 
time is thus indeed a factor ~ N larger than that for just 
two atoms. The phase of the atomic dipole diffuses slowly 
which is analogous to the slow phase diffusion of the field 
of a laser. This result is significant especially with an eye 
to possible applications as an ultrastable local oscillator 
or frequency reference. 

In summary, we have shown that steady state superra- 
diance can be achieved with atoms with an ultra-narrow 
transition in a cavity. Alkaline-earth-like atoms are prime 
candidates for realizing such systems and open up a new 
regime of CQED characterized by an extremely small 
critical photon number. In future work we will study 
the noise properties of the emitted light and the corre- 
lated atomic state. It will be intruiging to consider the 
crossover from steady-state superradiance to lasing that 
occurs as the system goes from the bad cavity limit to 
the good cavity limit. 
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